Majorana surface states of superfluid 3 He A- and B-phases in a slab 
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Motivated by experiments on the superfluid 3 He confined in a thin slab, we design a concrete 
experimental setup for observing the Majorana surface states. We solve the quasi-classical Eilen- 
berger equation, which is quantitatively reliable, to evaluate several quantities, such as local density 
of states (LDOS), mass current for the A-phase, and spin current for the B-phase. In connection 
with realistic slab samples, we consider the upper and lower surfaces and the side edges including 
the corners with several thicknesses. Consequently the influence on the Majorana zero modes from 
the spatial variation of /-vector for the A-phase in thick slabs and the energy splitting of the zero- 
energy quasi-particles for the B-phase confined in thin slabs are demonstrated. The corner of slabs 
in the B-phase is accompanied by the unique zero-energy LDOS of corner modes. On the basis 
of the quantitative calculation, we propose several feasible and verifiable experiments to check the 
existence of the Majorana surface states, such as the measurement of specific heat, edge current, 
and anisotropic spin susceptibility. 

PACS numbers: 67.30.ht, 74.45. +c 



I. INTRODUCTION 

Majorana quasi-particles (QPs) and Majorana 
fermions have attracted much attention in the wide 
research field, ranging from high energy physics to low 
temperature physics of ultracold atomic gases,-^ and for 
the application to topological quantum computations^ 
The Majorana QP and Majorana fermionic operator 



are defined by 7' 
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and \fd = iff, respectively, which 



imply that the particle and antiparticle are identical. 
It has been proposed that the Majorana nature brings 
new physics, such as non-abelian statistics of vortices 
in chiral superfluids^ and Ising-like spin dynamics for 
time-reversal invariant superfluidsJi — The Majorana 
nature itself is an intriguing subject to further studies. 

Candidate systems that exhibit the Majorana nature 
are quite rare, e.g., spin-triplet superconductors or su- 
pcrfluids, p-wave Feshbach resonated superfluids^'a frac- 
tional quantum Hall systems with the 5/2 filling in- 
terfaces between a topological insulator and an s-wave 
superconductor^ and s-wave supcrfluids with particu- 
lar spin-orbit interactions #ii 

Spin-triplet superconductors or supcrfluids satisfy the 
following Bogoliubov-de Gennes equation with the wave 
function tp n and the energy E n : 



f dr2U{ri,r2)<Pn(r 2 ) = E n (p n (n). 



(1) 



The eigenstates yield one-to-one mapping between the 
positive energy states ipE and the negative energy states 
ip -e = t x <Pe owing to the symmetry H = —t x H*t x , 
where t x is the Pauli matrix in the particle-hole space. 
The symmetry of the wave function leads to the relation 



„t 



Therefore, the 



of the Bogoliubov QP operator 7^ = 7i B . 
Bogoliubov QP with zero-energy is the Majorana QP, 
7q = 70 . Zero-energy bound states of the QPs appear 
whenever the underlying potential for the QPs changes 
its sign. The situation appears at edges, surfaces, or 



vortices with odd integer winding numbers, where the 
QPs have zero-energy Andrcev bound states. For the 
edge or surface, the Majorana fcrmion surface state = 
^> exists as well as 7g = 70^^ 

Among known possible spin-triplet superconduc- 
tors! 3 ^ UPt 3 , UGe 2 , and URhGe, Sr 2 Ru0 4 is a prime 
candidate of chiral spin-triplet superconductors. The 
precise pairing symmetry, however, has not yet been iden- 
tified and has been under strong discussion. A p x + ip y 
scenario is reexamined from various aspects^ Obviously, 
we want to experiment with more candidate materials for 
observing the Majorana nature. 

The superfluid 3 He consists of spin-triplet Cooper 
pairsJ^ For the 3 Hc there is a huge amount of exper- 
imental data by intensive researches for long yearsJ^ 
There are two stable phases, ABM- (A-) and BW- (B-) 
phases, for the superfluid 3 Hc in a bulk without a mag- 
netic field. The A-phase is stabilized in narrow regions at 
high temperatures and high pressures while the B-phase 
is stabilized in other wide regions within the superfluid 
phase. The order parameter (OP) of the A-phase has 
point nodes toward /-vector, which signifies the direction 
of the orbital angular momentum, or orbital chirality, 
whereas the B-phase has a full gap. The Majorana na- 
tures for the A- and B-phases differ due to the topology 
of the gap on the Fermi surfaced Time reversal sym- 
metry is also different for the A- and B-phases. Since 
up-up and down-down spin Cooper pairs have the same 
chirality for the A-phase, the time-reversal symmetry is 
broken and net mass current flows along the boundary or 
the edge. In contrast, since up-up and down-down spin 
Cooper pairs have the opposite chirality, the B-phase is 
a time-reversal invariant. Instead of the mass current 
canceled by up-up and down-down spin Cooper pairs, 
the spin current flows. The difference of time reversal 
symmetry between the A- and B-phases is analogous to 
the difference between the quantum Hall state and the 
quantum spin Hall stated 
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Our aim is to propose a concrete experimental design 
to detect the Majorana nature based on quantitative cal- 
culations and understand systematically how to observe 
the Majorana nature in the superfluid 3 He A- and B- 
phases. The quantitative calculations are performed by 
quasi-classical theory which yields information of QPs. 
The quasi-classical theory is valid when A/Ep <C 1, 
where A is superfluid gap and Ep is the Fermi energy. 
Since A/Ep ~ 10~ 3 , the theory is appropriate for the su- 
perfluid 3 Hei22 The quasi-classical framework is well es- 
tablished within the weak-coupling limit. The superfluid 
3 He at the lower pressure can be described within the 
weak-coupling scheme^ without delicate strong-coupling 
corrections. 

In the slab geometry, both the A- and B-phases are 
stabilized by changing thickness and temperature even at 
zero pressure i 24 ' 25 Therefore, we can deal with the A- and 
B-phases by the quasi-classical theory within the weak- 
coupling limit. Here, we introduce several systems of the 
slab geometry. By Bennett et ai.;2& the superfluid 3 Hc is 
confined in a thin slab box with a thickness D = 0.6 /xm 
and an area of the base 10 mm x 7 mm. The thickness is 
of the order of 10£o at P = 0, where £o is the coherence 
length at T = 0. At ISSP, University of Tokyo^ sample 
disks with a thickness D = 12.5 /im, which is of the 
order of the dipole coherence length, and a diameter 3 
mm are used. The ISSP group can rotate the sample 
disks with the highest speed in the world and investigates 
half-quantum vortices. At RIKEN^ the inter-digitated 
capacitors are used for making a film of the superfluid 
3 He. The group can control thicknesses of the film from 
0.3 to 4 /im. 

Under those experimental situations, we consider the 
surface states from not only the upper and lower surfaces 
in a slab but also the side edges including the corners 
and evaluate the dependence of the Majorana state on 
the thickness of a slab. By the quasi-classical calcula- 
tion, we show the"Dirac valley" (Fig. [TJa)) for the A- 
phase at the side edges. In a thick slab, we consider the 
fc z -component of the OP, which was neglected in the pre- 
vious work 2 ^ since Z-vector was assumed to be perpendic- 
ular to the thin slab. Several groups have demonstrated 
that the Majorana fermion surface state has the disper- 
sion of E = (A/kp) |fc|| | with the surface perpendicular 
to the z-axis in the B-phase, where kp is the Fermi wave 
number and fen = (k x ,k y )£-~— This implies that the sur- 
face state consists of a single Majorana cone. We show 
the Majorana cone (Fig.QTb)) by the quasi-classical cal- 
culation for a thick slab and investigate the variation of 
the dispersion by the thickness of a slab (Fig. [TJc)). The 
zero-energy state at the corner is also discussed. 

The arrangement of this paper is as follows: In Sec.lTfl 
we formulate the quasi-classical theory based on the 
quasi-classical Green's function which gives quantitative 
information of QPs. We explain numerical methods in 
Sec. IIII1 which is supplemented by the symmetry consid- 
erations of the quasi-classical Green's function given in 
Appendix. In Sees . II VI and fVl for the A- and B-phases, we 




FIG. 1. (Color online) Stereographic views of typical disper- 
sion relation, (a) Dirac valley for the A-phase in a thin slab 
D = 8£o at the left side edge "L" and the right side edge "R" . 
(b) Majorana cone for the B-phase in a thick slab D = 30£o 
at the upper or lower surface, (c) Split Majorana cone for 
the B-phase in a thin slab D — 14£o at the upper or lower 
surface. The units of energy and wave length are ttUbTc and 
kp, respectively. 



show results of the spatial structures of the OP, the cur- 
rent distribution, and the local density of states (LDOS) 
for QPs which relates to the Majorana modes. In Sec. lVIl 
we discuss difference of the Majorana zero modes in the 
A- and B-phase and propose several experimental designs 
to observe the Majorana nature. In addition, we discuss 
the Majorana zero mode in the stripe phase. The final 
section is devoted to a summary. 



II. QUASI-CLASSICAL THEORY AND ORDER 
PARAMETERS 

We start with the quasi-classical spinful Eilenberger 
equation,— which has been used for studies of the super- 
fluid 3 He i 25 ' 31 i 32 The low energy excitation modes at the 
surface are discretized in the order of A/(kpL)^ where 
L is a length of the system. If L is a macroscopic length, 
which is much larger than fc^ 1 , the low energy excita- 
tion modes at the surface can be regarded as continuous 
spectra. Because there is the exact zero-energy excitation 
at the surface of the superfluid 3 Hc A-phas e 12 ' 29 and B- 
phase^— the quasi-classical theory can be used to discuss 
the Majorana QP. 

The quasi-classical Green's function g(k,r,uj n ) is cal- 
culated using the Eilenberger equation 



ihv(k) ■ Vg(k,r,u n ) 



A(fc,r)t 



(2) 



where the "ordinary hat" indicates the 2x2 matrix in 
spin space and the "wide hat" indicates the 4x4 ma- 
trix in particle-hole and spin spaces. The quasi-classical 
Green's function is described in particle-hole space by 



g(k,r,u) n ) 



g(k,r,uj n ) if(k,r,u> n )\ 
-if(k,r,w n ) -g(k,r,u n )J ' 



with the direction of the relative momentum of a Cooper 
pair fe, the center-of-mass coordinate of a Cooper pair r, 
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and the Matsubara frequency u> n — (2n + l)7rfcsT. The 
quasi-classical Green's function satisfies a normalization 
condition g 2 = — ir 2 l. The Fermi velocity is given as 
v(k) = vpk on the three dimensional Fermi sphere. 

We solve Eq. © by the Riccati method^2M4 We 
introduce Riccati amplitude a = (1 + 9) f and b = 
(1 + g)^ 1 f related to particle- and hole-like projections 
of the off-diagonal propagators, respectively. Equation 
([2|) can be rewritten as Riccati equations 

hv(k) ■ Vo(fe, r,u„) 

= A(fc,r) - a(k, r, w„)A(fc, r) f a(fc, r, w„) 

- 2w n a(fc,r, w n ), 



— hv(k) ■ V£>(fe, r,a; n ) 

= A(fc, r)t - S(fc, r, w n )A(k, r)b{k, r, u n ) 

-2w n 6(fe,r,a; n ). (4) 

The equations are solved by integration toward fc for 
a(k,r,u) n ) and toward — fe for b(k,r,u! n ). From the 
Riccati amplitude, the quasi-classical Green's function 
is given as 



(l + aby 






(i + Sa)- 



l-ab 



2ia 



-2ib -(1-ba) 



(5) 



The self-consistent condition for the pair potential 
A(fc, r) is given as 



A(fc,r) = N 7rk B T (V(k,k')f(k',r,u r , 

— U} c <UJ n <.0J C 



k' 

(6) 



where Nq is the density of states in the normal state, 
lu c is a cutoff energy setting w c = AOirksTc with the 
transition temperature T c in a bulk, and (•••)& indi- 
cates the Fermi surface average. The pairing interaction 
V(k, k) = 3gik ■ k' for Cooper pairs with an orbital an- 
gular momentum L = 1, where g\ is a coupling constant. 
In our calculation, we use a relation 



1 =] T 



27rfc B T V — 



(7) 



Spin-triplet OP is defined by a vectorial notation 

A(k,r) = (i&a y ) ■ A(fc,r), (8) 

with the Pauli matrix <r in spin space. The complex 
vector A is perpendicular to the spin S of a Cooper pair, 
namely, A • S = 0. The A can be expanded in orbital 
momentum directions, 



A Al (fc,r) = Apifr)}* 



(9) 



where ^^(r) is a complex 3x3 matrix with a spin index 
/i and an orbital index i. The repeated index implies 
summation over x, y, and z. 

This paper is discussed for the chiral state in the A- 
phase, polar state, B-phase, and planar state. Symbolic 
descriptions of the OP in their state are the following: 
For the chiral state, 



for the polar state, 



for the B-phase, 



dz \k x iky 



A — d z k Xl 



A — d x k x + dyky + d z k z , 



and for the planar state, 

A = d x k x dyky, 



(10) 

(11) 
(12) 

(13) 



where d- vector is perpendicular to the spin of a Cooper 
pair, ci-vector and the momentum direction are per- 
mitted to rotate spherically in the spin and momentum 
space, respectively. 

By using the self-consistent quasi-classical Green's 
function, the mass and spin currents are calculated by 

j(r) = mN nk B T ^2 ( v ( k ) Im [9o(k, r, w„)])fe, 

(14) 



-LO c <~~UJ n <~~CJ c 



3^(r) = ^N TTk B T Y, (v(k)lm[ 9tt (k,r,uj n )}) k , 

— LU c <UJ rl <UJ C 

(15) 



respectively, where m is the mass of the 3 He atom and 
g^ t is a component of the quasi-classical Green's function 
g in spin space, namely, 
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.90 + 9z 9x - i9 y 



J)x + Wy 90 - gz 

LDOS for energy E is given by 
N(r,E) = (N(r,E,k)) k 

])*. (16) 

where 77 is a positive infinitesimal constant and 
N(r, E 7 k) is angle-resolved LDOS. Typically, we choose 
r) = 0.0037rfcsT c . For obtaining go(k,r ,u n )\iu n ->E+iri, 
we solve Eqs. (U) with r) — iE instead of uj n under the 
pair potential obtained sclf-consistently. 



III. SYSTEM GEOMETRY AND NUMERICAL 
METHODS 

We consider a cross-section A of a slab with a thickness 
of D along the ^-direction and a macroscopic length of L 
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FIG. 2. (Color online) (a) A schematic configuration of a slab 
with a thickness of D along the z-direction and a macroscopic 
length of L along the x-direction. We consider a cross-section 
A (0 < x < L,0 < z < D). (b) Indexes of positions in 
the cross-section A. We label {x, z) = (0, 0), (0, D/2), (0, D), 
(L/2,0), and (L/2, D/2) as "1", "2", "3", "4", and "5". Paths 
z = 0, z = D/2, z = D, and z = 0.5£o along the z-axis are 
labeled as "I", "II", "III", and "IV". 



along the x-direction, shown in Fig. HJa) . It is assumed 
that the quasi-classical Green's function and OP are ho- 
mogeneous along the y-direction. We also assume that 
the surfaces of the slab are specular, where specularity is 
controlled by coating the surface with 4 He atoms £2- Un- 
der the boundary condition, the quasi-classical Green's 
function or the Riccati amplitude changes only the di- 
rection of the relative momentum by mirror reflection at 
a surface H surf , namely, g(k, R surt , u n ) = <?(fc, R sm i, w«) 
or a(k,R suri , ui n ) = a(fc, -R surf , w„) and 6(fc, R suli , u) n ) = 
b(k, Rs UT {, w n ), where k = k — 2n(n-k) with a unit vector 
n which is perpendicular to the surface. 

We solve Eqs. (Q| by the numerical integration to- 
ward the fc-direction for a and toward the —fc-direction 
for b from the surface to the position r. Since 
a(fc, -R S urf, w„) = a(k, R sul -[ , uj n ) , the initial value of a 
at a surface i? sur f is determined by the numerical inte- 
gration toward the fc-direction from the other surface to 
-Rsurf. Eventually, we have to solve the numerical inte- 
gration along the sufficiently long path by changing the 
direction from k to fc at the surface so that the integra- 
tion path gives the same value of a under arbitrary ini- 
tial values, and also b. Since the integration path changes 
the direction at the surface, simple calculation is difficult. 
For overcoming the difficulty, we exploit that the change 
in the direction of the integration path from fc to fc can 
be regarded as the change in the relative momentum of 
the pair potential from k to fc. For example, we can 
substitute —An X for An X in the pair potential instead of 
the reflection of the integration paths at a surface x = 0. 
Similarly for the other surfaces, we can integrate Eqs. (|4]) 
along the straight long paths in the region connected to 
the antiperiodic pair potential infinitely. 

We calculate the quasi-classical Green's function in a 
range -L/2 < x < L/2 and -D/2 < z < D/2 instead 
of a cross-section A (0 < x < L, < z < D) to reduce 
computational time. We show the reduction method in 
Appendix in addition to other reduction methods with 
the symmetry of the quasi-classical Green's function. 



IV. A-PHASE 

We discuss here the situation in which the thickness of 
a slab D is changed, where the A- and B-phases are sta- 
ble in thin and thick slabs, respectively^ We present the 
calculated results of OP, LDOS, and mass current for the 
A-phase and spin current for the B-phase at T — 0.2T C . 
In each phase, the results depend on the thickness D. 
The length of the cross-section A is taken as L = 40£o 
throughout the paper where the OP is recovered to the 
bulk value at the center of the system. Note that the 
following results are unchanged for the length L longer 
than 40£o- In addition, we consider the low pressure 
limit P — > within the weak-coupling limit. The spa- 
tial positions of the results are indicated by the indexes 
in Fig. Hlb). In this section and next section, we use 
the units £o = Flvf /2i:kBT c , irkBT c , and Nq for length, 
energy, and LDOS, respectively. 

The bulk of the A-phase is in the chiral state. The OP 
in the chiral state is described by-i2 

= d^(m + in)i, 

where m-vector and n-vector are perpendicular to each 
other. Z-vector, which signifies the orbital chirality, is de- 
fined as I = m x n. At the edge of a slab, since the 
normal orbital component to the edge vanishes, the po- 
lar state will be realized. The OP in the polar state is 
described by-i£ 

Afii = dftiTii, 

where we regard n-vector as perpendicular to the edge. 

In a slab cell used by Bennett et al££- where the thick- 
ness is 0.6 /im w 8£o, it is expected that the A-phasc 
is stable at T = 0.2T C £H Since the thickness D is much 
shorter than the dipole coherence length ~ 1000£o, d- 
vector, which characterizes the spin state of the OP in 
the A-phase, points to the z-direction in the absence of a 
magnetic field as long as Z-vector is parallel to the z-axis. 
Z-vector is parallel to the z-axis everywhere except near 
the side edges at x = and x = L in the slab shown 
in Fig. [5] The length scale of the spatial variation of l- 
vector is short, whose order is the coherence length. On 
the other hand, d-vector can not vary spatially since the 
order of the length scale is the dipole coherence length, d- 
vector is spatially uniform even under a magnetic fields 
In this paper, we fix the direction of d-vector to the z- 
axis. 



A. D = 8£o : Thin slab for A-phase 

For D = 8, the OP is described by A z (k,x) = 
A zx (x)k x + A zy (x)k y , where the relative phase between 
A zx and A zy is ir/2. The OP is uniform along the z- 
direction and varies along only the x-direction. Since 
the slab is thin, the A: z -component of the OP is sup- 
pressed. Note that the uniformncss along the z-direction 
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and the suppression of the fc z -component are not sup- 
positions like in the previous work^ but results by the 
self-consistent calculation. The profile of the OP along 
the x-axis is shown in Fig. Eta). Because of the spec- 
ular boundary condition, the fc x -component perpendicu- 
lar to the edge becomes zero at x = and x = L. In 
contrast, the parallel fcj,-component is enhanced by com- 
pensating for the loss of the fc^-component at the edge, 
where the polar state is realized. The fc^-component 
increases and the fc y -component decreases toward the 
bulk region around x = L/2; thus, the chiral state 
with k x + iky is attained. We can construct Z-vector as 
li = —ieiji;A* Z jA z k/\A\ 2 , where e^fc is the totally anti- 
symmetric tensor and |A| 2 = A* zi A Z i is the squared am- 
plitude of the OP. Z-vector points to the z-direction in 
the bulk region and vanishes at the edge. Since the chi- 
ral state is realized except at the edge, the Majorana 
fermion edge state = \f r exists at side edgesi&12 

The mass current j y (x) is shown in Fig. EJb). The 
mass current flows circularly along the side edge of the 
slab. Experimental values^ are used as coefficients of 
Eq. (|14[) so that a quantitative value of the mass current 
is obtained. By applying a magnetic field perpendicular 
to a slab, we can produce a spin imbalance due to the 
Zeeman shift between the up-up spin pairs and the down- 
down spin pairs. This results in a net spin current in 
addition to the mass current. 

Figure [3£c) shows LDOS at the edge x = (the po- 
sition "2" in Fig. [2]) and the bulk x = L/2 ("5"). It 
is clearly seen from the line "2" for x = that the 
LDOS with a substantial weight appears at a zero-energy, 
corresponding to the Majorana edge mode, because the 
distance L between the edges is macroscopic.^^ This 
implies that the LDOS is expressed as N(E, x = 0) = 
N(E = 0,x = 0) + aE 2 in the vicinity of E = 0. The 
first term comes from the Majorana edge mode, and the 
second term comes from the point node of the chiral state 
in the bulk A-phase. The peak at E ss 0.65 comes from 
the gap of the chiral state in the bulk and the peak 
at E ss 0.8 comes from the gap of the polar state at 
the edge. The LDOS at x = L/2 (the line "5") shows 
a typical behavior of the point node spectrum, namely, 
N(E, x = L/2) oc E 2 . In Fig. G^d) we show the extent 
of the Majorana edge mode at E = toward the bulk 
from the edge at x = 0, which spreads over the order of 
5£o- The spectrum of the edge gradually changes into the 
bulk spectrum. 

The LDOS N(r, E) is obtained by averaging the angle- 
resolved LDOS N(r,E,k) on the Fermi surface, which 
has the peak at the energy of the surface Andreev bound 
state. The dispersion relation of the surface Andreev 
bound state can be evaluated by following the peak in 
each momentum direction. The angle-resolved LDOS 
N(E, x = Q,6) with = 0° and N(E, x = 0, 0) with 8 = 
90° are shown in Figs.|3te) and[3jf), respectively, where 8 
is the polar angle from the fc z -axis and <f> is the azimuthal 
angle from the fc^-axis on the three-dimensional Fermi 
sphere. The angles and momentum directions have the 



relations: tan 8 = ^Jk 2 + k 2 /k z and tan</> = k y /k x . The 

angle-resolved LDOS is symmetric between 0° < 8 < 90° 
and 180° > 8 > 90°, and between -90° < cj) < 90° and 
270° > <fi > 90°. The zero-energy peaks of the angle- 
resolved LDOS appear if </> = 0° except for 8 = 0°. 
The zero-energy peaks come from the zero-energy An- 
dreev bound state of QPs at the edge. The QPs without 
the fc^-component have the zero-energy Andreev bound 
state because they feel the 7r-phase shift of the pair po- 
tential for the relative momentum from (k x ,k y = 0,k z ) 
to (—k x ,k y = 0, k z ). This is clear from the Eilenberger 
equation @ since the quasi-classical Green's function 
g(k) depends on only the fc-component of the pair po- 
tential k x +ik y . The angle- resolved LDOS is asymmetric 
with respect to E = (Fig. CJf)) because the supcrfluid 
state in the A-phase has the chirality^S The asymmetric 
angle-resolved LDOS is related to the mass current. 

The peak energy of angle-resolved LDOS at the side 
edges as a function of k z for k y = and as a function of 
k y for k z = is shown in Figs.EHg) and|3£h), respectively. 
It is found that the energy is dispersionless along the k z - 
axis and two linear dispersions appear near E = along 
the k y -axis, where each dispersion comes from the left 
(x = 0) and right (x = L) edges. Since the directions 
of the mass current are different both edges, the linear 
dispersions have the opposite slope. The linear Dirac 
dispersions continue along the fc z -direction, which form a 
"Dirac valley" . The results are consistent with the strict 
quantum level structures of the Majorana fermion edge 
state by Bogoliubov-de Gennes equation^ 

B. D = 14£ : Thick slab for A-phase 

In a thick slab, the /c z -component of the OP is induced. 
For D = 14, the fc z -component appears near the edges at 
x = and x = L to avoid the polar state except at z = 
and z = D, shown in Figs.2ta)-U£c). The OP is described 
by A z (k,r) = A zx (r)k x + A zy (r)k y + A zz (r)k z , where 
the phases of A zz and A zx are the same and the relative 
phase between ^4 ZZ and A zy is ir/2. Since the thickness of 
D = 14 is small, the A zz does not recover the bulk value 
of A zx at the side edges x = and x = L along z = D/2 
(Fig. IDJc)). As similarly at the edges z = and z = D, 
Z-vector is perpendicular to the edges at x = and x = L 
by the induced fc z -component, shown in Fig. |U[d). 

The mass current j y (r) is shown in Fig. [5ta), and its 
profiles along "I" (z = 0), "II" (z = D/2), and "III" 
(z = D) are shown in Fig. [Sfb). The mass current flows 
clockwise relative the local Z-vector direction. The j y 
from Z-vector in the bulk region flows to the negative 
(positive) y-direction at x = (x = L). At the side 
edges, since Z-vector points to the negative a;-direction, 
the j y from the edges flows to the negative and positive y- 
directions in z < D/2 and z > D/2, respectively. There- 
fore, the edge mass current is strengthened in z < D/2 
(z > D/2) and weakened in z > D/2 (z < D/2) at the 
edge x = (x = L). The symmetry of the mass current 
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FIG. 3. (Color online) Calculated results for the A-phase with D = 8. Profiles of OP (a), and mass current j y (b) along the 
x-axis. (c) LDOS N(E,x) at "2", and "5" defined in Fig.^b). (d) Zero-energy LDOS N(E = 0,x) from the edge at x = to 
the bulk at x = L/2. Angle-resolved LDOS N(E,x = 0,9) with <j> = 0° (e), and N(E,x = 0,0) with 9 = 90° (f). The peak 
energy of angle-resolved LDOS at the side edges as a function of k z for k y = (g), and as a function of k y for k z = (h), where 
"L" ("R") is for the left (right) edge. In this and the following figures, the units of length, energy and A^i, LDOS, and wave 
length are £o, ^ksT c , No, and fc_F, respectively. 




x 



FIG. 4. (Color online) The OP for the A-phase with D=14. 
(a) Amplitude of each component. Profiles of each component 
along "I" (b), and along "II" (c). (d) Orientations of I- vector 
texture in the x-z plane of a slab. 



between z < D/2 and z > D/2 is broken despite the 
symmetry of the OP. 

Figure EJc) shows LDOS at (x, z) = (0,0), (0,D/2), 
(0,.D), and (L/2, D/2). At (0,0) and (0, D) (the lines 
"1" and "3"), there are the zero-energy LDOS and a 
small peak at a low energy indicated by the arrow in 
Fig.JSfc). Although L vector points to the edge, the zero- 
energy LDOS also exists at (0,L>/2) (the line "2"). The 



lines "1" and "3" are slightly decreased and increased, 
respectively, in low energy from the small peak as energy 
increases. The line "2" is almost constant in low energy. 
The LDOS (the lines "1", "2", and "3") are enhanced 
sharply near E « 0.5. The LDOS at (L/2, D/2) (the 
line "5" ) shows the typical spectral behavior of the point 
node. In Fig. [5jd) we show the extent of the zero-energy 
edge mode toward the bulk from the edge at x = 0. 
The amounts of the zero-energy LDOS at x = along 
z = D/2 and z = D (the lines "II" and "III") arc smaller 
than that along z = (the line "I"). The extent of the 
zero-energy LDOS along "III" is of the order of 5£o as 
in the thin slab D = 8; however, the zero-energy LDOS 
along "I" and "II" are more widely extended. 

The angle-resolved LDOS N(E,r,6) with <j) = 0° at 
(0,0), (0, D/2), and (0,D) are, respectively, displayed 
in Figs. (He), EJf), and[5Jg). The zero-energy peaks of 
the angle-resolved LDOS appear except for 8 = 0°. The 
peaks of the line 9 = 0° at the finite low energy ap- 
pear from the upper and lower surfaces because the k z - 
componcnt of the pair potential is non- vanishing. Since 
thickness D = 14 is small, the original zero-energy peak 
of 9 = 0° perpendicular to the upper and lower surfaces is 
split into the finite energy peaks, which is discussed par- 
ticularly for the B-phasc in next section. The splitting 
peaks at (0, 0) are shown in Fig.^h) only at k z = ±1 for 
k y = 0. The finite energy peaks compose the small peak 
of the LDOS in the lines "1" and "3" in Fig.^c). Since 
the position (0, D/2) is apart form the upper and lower 
surfaces, the energy peaks of 9 = 0° are small. Conse- 
quently, the peak of the LDOS at low energy is not clear 
(the line "2" in Fig. [5(c)). The point nodes of the OP 
near the side edges have a tilt to the x-direction accord- 
ing to ^-vector, or the anti-nodes of the OP lie to some 
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angles from the £; z -axis. Thus, the incident QPs from 
(0, 0) with low angles from the fc 2 -axis do not feel the 
clear small gap in the vicinity of the point nodes. For in- 
stance, the line 9 = 30° in Fig. (SJe) should be compared 
with that in Fig. (3je) for the thin slab D = 8. There- 
fore, the gap-like LDOS enhanced sharply near E ps 0.5 
appears (the line "1" in Fig. E[c)). On the other hand, 
the incident QPs from (0,D) with low angles from the 
fcz-axis are reflected at the edge and regarded as with 
low angles from the — fc 2 -axis. Since the QPs feel the 
small gap in the vicinity of the point nodes, they have 
the clear peaks in the angle-resolved LDOS (Fig. EJg)). 
Therefore, the LDOS is slightly increased in low energy 
as energy increases (the line "3" in Fig. (5jc)) . Since the 
QPs at (0, D/2) have the two kinds of the characteristic 
momentum, the LDOS is almost constant in low energy 
(the line "2" in Fig. Etc)). 

The LDOS at the side edge is different and depends 
on the thickness of the slab. In the thin slab (D = 8), 
l-vector points to the z-direction everywhere and the 
LDOS is N(E) = JV(0) + aE 2 in low energy. In con- 
trast, in the thick slab (D = 14), I- vector points to 
the x-direction near the side edge and the LDOS is 
N(E) = N(0) + N gap (E), where N gap (E) is almost zero 
in low energy and increases sharply near E sa 0.5 like the 
LDOS with the full gap state. The edge mass current 
varies along the side edge in the thick slab, which is also 
a great difference with the thin slab. 

Note that the A-phase is metastable in the slab with 
D = 14 at T = 0.2T C in the weak-coupling limit for 
P 0M- In the same slab, the A-phase is stabilized 
at T = 0.9T C )2£ however, Z-vector points to only the z- 
direction, because the fc z -componcnt is absent owing to 
the longer coherence length at higher temperature. The 
mechanism is similar to that the A-phase is more sta- 
ble than the B-phase, which has the k x -, k y -, and k z - 
components, in a thin slab at low temperature. The cal- 
culated results shown in this subsection is not in free 
energy minimum at P — 0; however, the strong-coupling 
effect by pressure stabilizes the A-phase texture. Thus, 
the above characteristics from the texture will be ob- 
served at finite pressure. 



V. B-PHASE 

The B-phase is stable when the thickness of a slab is 
larger than w 13£ at T = 0.2T C ^ The OP of the B- 
phase is described byi£ 

(A xx \ 
A^i = R^(n d ,9 d ) A yy , 
V A ZZ J 

where Rp,i(rid,9d) is a rotation matrix with a rotation 
axis rid an d a rotation angle 0d about rid- The rota- 
tion matrix gives the relative angle between the orbital 
momentum and the direction along which the spin of a 



Cooper pair is zero. The spin state is stable by the dipolc- 
dipole interaction when 9 d = 9 L = cos _1 (— 1/4) and rid is 
perpendicular to the surface in the absence of a magnetic 
fields If the thickness of a slab is much smaller than the 
dipole coherence length, n d is locked to the z-axis. Thus, 
we derive A xx (r), A yy (r), and A zz (r) with the uniform 
rotation matrix R{z,9l). Note that the three compo- 
nents have the same phase. 

At the surface of a slab, since the normal component to 
the surface of the orbital state vanishes, the planar state 
will be realized. The OP in the planar state is described 
bylS 

(A xx 0\ 
A„i = R„i(n d ,9 d ) A yy , 
\ 0/ 

where we regard the z-axis as perpendicular to the edge. 
If 9 d = 0, the OP can be written by other descriptions as 

where A n = -(A x - iA y )/V2, A u = (A x + iA y )/V2, 
and A p = A xx = A yy . Therefore, in the planar state, 
S z = +1 spin state has L z = — 1 orbital angular momen- 
tum and S z = — 1 spin state has L z = +1 orbital angular 
momentum. 



A. D = 30£o: Thick slab for B-phase 

We consider that the thickness of a slab is much longer 
than the coherence length, namely, the thickness can be 
regarded as macroscopic for the OP, and much shorter 
than the dipole coherence length. The requirement is 
satisfied for a thickness of D = 30. It is clearly seen 
that the component of the OP perpendicular to the edge 
becomes zero in Fig. Eta). The polar state k y occurs at 
the corner "1" [x = 0, z = 0). The planar state k x — ik y 
for the up-up spin Cooper pairs and k x + ik y for the 
down-down spin Cooper pairs is attained in the middle 
region around "4" (x = L/2, z = 0), shown in Fig. E^b). 
The planar state k y ± ik z is realized at the side edge "2" 
(x = 0, z = D/2) and the OP is recovered to the bulk 
value of the B-phase around "5" (x = L/2,z — D/2), 
shown in Fig. EJc). At all edges (x = 0, x = L, z = 0, 
and z = D), the Majorana fermion edge state exists^— 

The spin current j^ y {v), jl y (v), (r), and j v sz {r) are 
shown in Fig. [TJa). Since the spin current flows in the 
three-dimensional real space, we show the schematic flow 
of the spin current in Fig. [T^b). The spin current for the 
i-componcnt of spin turns around the i-axis. This is un- 
derstandable as follows: Since the up-up spin Cooper 
pairs have negative chirality and the down-down spin 
Cooper pairs have positive chirality, the z-component of 
the spin current turns around the z-axis. Other compo- 
nents are also understood in the same manner. Note that 
the schematic flow in Fig. E^b) is derived from An before 



8 



(a) 



j y [kg m- 2 s- 1 ] 








i5 


♦ 1 


3^ 
















60° 


90° 


4ft 









40 " 
(g) K 



0.4 ^ M 





E 



(i 




FIG. 5. (Color online) Calculated results for the A-phase with D=14. (a) Mass current j y (r) in the x-z plane, (b) Profiles of 
the mass current j y (r) along paths "I", "II", and "III", (c) LDOS N(E,r) at "I", "2", "3", and "5". (d) Zero-energy LDOS 
N(E = 0,r) along "I", "II", and "III". Angle-resolved LDOS N(E,r,8) with cj> = 0° at "1" (e), "2" (f), and "3" (g). (h) The 
peak energy of angle- resolved LDOS as a function of k z for k y = at "1" . 




FIG. 6. (Color online) The OP for the B-phase with D=30. 
(a) Amplitude of each component. Profiles of each component 
along "I" (b), and along "II" (c). 




rotation by R(z,8l). The rotation axes for the x- and 
y-component of the spin current are changed correspond- 
ingly. 

Figure (5Ja) shows the LDOS. It is seen that only the 
line "1" (x = 0, z = 0) has the zero-energy LDOS from 
the zero-energy Andreev bound state of QPs without the 
^-component. The lines "2" and "4" (x = 0, z = D/2 
and x = L/2,z = 0) are almost the same and show the 
linear relation N(E, r) oc E near E = which reflects 
the surface Andreev bound states! 38 ' 39 The line "5" (x = 
L/2, z — D/2) corresponds to the bulk LDOS where the 
full gap is expected. In Fig. [SJb) we show the extension 



FIG. 7. (Color online) (a) The spin current for the B-phase 
with D — 30. The unit of the spin current is 10 -9 J/m 2 . 
(b) A schematic flow of the spin current for S x , S y , and S z 
components. 



of the zero-energy state toward the middle from the edge 
at x = 0. The zero-energy LDOS of the line "I" (z = 0) 
reduces sharply as approaching the middle region from 
the edge. It is found that the value N(E = 0, x = 0) 
in the line "IV" (z — 0.5) is much smaller than that 
of the line "I". The line "II" along z = D/2 shows the 
absence of the zero-energy LDOS. Therefore, we conclude 



9 




hi ki Ki 



FIG. 8. (Color online) The LDOS for the B-phase with D = 
30. (a) LDOS N(E,r) at "1", "2", "4", and "5". (b) Zero- 
energy LDOS N(E = 0,r) along "I", "II", and "IV" (z = 
0.5). The peak energy of angle-resolved LDOS as a function 
of k y for k z = and as a function of k z for k y = at "2" (c), 
and as a function of k x for k y — and as a function of fc H for 
k x — at "4" (d). (e) and (f) are the same plots as (c) and 
(d), but at "1". 



that the zero-energy LDOS is localized at the corner of 
the order of £o- The peak energy of the angle- resolved 
LDOS is shown in Figs. [8jc) for the side edge and [8(d) 
for the lower surface. Figure [5(c) is as a function of k y 
for k z — and as a function of k z for k y = at "2" 
(x = 0, z = D/2). Figure [5(d) is as a function of k x 
for fcj, = and as a function of k y for fc x = at "4" 
(x = L/2,z = 0). These show that the side edge and 
lower surface have the same dispersion relation for low 
energy. The linear dispersion forms the Majorana cone 
at the surfaced 

The peak energy of the anglc-rcsolved LDOS at the 
corner (0,0) is shown in Figs. [5(c) and [5(f). The energy 
is dispcrsionlcss along the k z - and fc^-axes for k y = 0. 
The Dirac valley-type dispersion is composed of the zero- 
energy LDOS at the corner (the line "1" in Fig. [5Ja)). 
This is different from the Majorana cone at the surface. 



B. D = 14£ : Thin slab for B-phase 

We consider D = 14 in thickness of a slab which is the 
critical thickness for stabilizing the B-phase below which 
the stripe phase becomes stabled Since the thickness is 
short, the fc z -componcnt of the OP can not completely 
recover the bulk value. Except that the fc z -component is 



FIG. 9. (Color online) The same as in Fig. [8] but for the 
B-phase with D = 14. 



reduced, the qualitative features of the OP and the spin 
current are the same in the case of D = 30; however, the 
LDOS shows different spectrum. 



Figure [9] shows the LDOS and the peak energy of the 
angle- resolved LDOS, which are compared with those in 
Fig.[8j The differences between them are following: Since 
the thickness of the slab is short, the line "5" for the 
center (x = L/2,z = D/2) has the spectrum from the 
full gap and the linear relation near E = 0, shown in 
Fig. [DJa). The latter comes from the contributions ex- 
tending from the upper and lower surfaces. The line "2" 
for the side edge (x = 0, z = D/2) in Fig. |H(a) has the 
small peak at E ps 0.35 composed of the deformed Ma- 
jorana cone in 0.5 J$ \k z \ ^ 0.7, shown in Fig. H(c), be- 
cause the QPs with the /c z -component reflect that the 
A: z -component of the OP does not recover the bulk value. 
The peak energy at "4" for the lower surface in Fig. [9(d) 
implies that the lowest energy of the Majorana cone is 
lifted from zero because of the tunneling of the zero- 
energy modes bound at two surfaces corresponding to 
the upper and lower surfaces. For the same reason, the 
zero-energy modes at the corner is split for k z = ±1 in 
Fig. [9(e) and for k y = or k x = in Fig. [9(f). By 
the energy splitting, the zero-energy LDOS is slightly re- 
duced; however, the extension of the zero-energy state is 
the same for the thick slab D = 30, shown in Fig. [5(b). 
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FIG. 10. (Color online) Calculated results for the B-phase at 
"4". (a) The peak energy of angle-resolved LDOS as a func- 
tion of k x for k y = and as a function of k y for k x = with 
several thicknesses D = 30, 12, and 8. (b) The dependence 
of the energy splitting on the thickness. The dashed line is 
proportional to exp( — D/3£o). 



Dependence of the energy splitting on the 
thickness 



The lowest energy of the Majorana cone at the lower 
surface is more lifted from the zero-energy as the thick- 
ness is shorter. If the thickness is macroscopically long, 
the zero-energy modes appear, shown in the left panel of 
Fig. flUTa) for the B-phase at "4" (x = L/2, z = 0) with 
D = 30. When the thickness is short, the zero-energy 
modes bound at two surfaces are hybridized with each 
other. They form the symmetric and anti-symmetric 
states with the opposite sign energy E + and re- 
spectively, where \E + \ = |£L| = £/ S pHt- The repre- 
sentative dispersion relation with D = 12 is shown in 
the middle panel in Fig. [TUta). In the case of more 
shorter thickness, the B-phase becomes the planar phase, 
which is metastable against the A-phase, because the k z - 
componcnt of the OP vanishes. Since the Fermi surface 
has the point nodes toward the direction of k x = k y = 0, 
the angle-resolved LDOS toward the direction corre- 
sponds to that of the normal state without the peak en- 
ergy. We show the dispersion relation with D = 8 except 



the point k x = k y = in the right panel in Fig HOf a) . 
There is a difference between the dispersion relations for 
D = 12 and D = 8 also in long wavelength which are 
curved and linear, respectively. 

The dependence of the energy splitting on the thick- 
ness is shown in Fig. [TUTb). The energy splitting has the 
relation S 8 plit ^ exp(— D/3£o) ; where 3£o is effective co- 
herence length. The exponential suppression is similar to 
the energy splitting of the Majorana zero-energy modes 
bound at two vortices in a chiral p-wave supcrfluid i 40 ! 41 



VI. DISCUSSION 



A. Majorana zero modes in the A- and B-phases 



In the A-phase which has point nodes, since the zero- 
energy QPs appear on the line k y — in the momen- 
tum space, the zero-energy QPs are dispcrsionless along 
the fc z -direction (Fig. [2Ig))- Therefore, the dispersion 
of the QPs forms the Dirac valley and the zero-energy 
LDOS spreads over the order of 5£o from the side edge 
(Figs. [3Jc) and[3][d)). In contrast, since the OP in the 
B-phase has full gap on the Fermi surface, only the QPs 
with the momentum perpendicular to the edge are the 
zero-energy Majorana QPs. Therefore, the dispersion of 
the QPs forms the Majorana cone (Figs. [5|c) and^d)). 
However, at the corner, the dispersion of the QPs forms 
the Dirac valley and the zero-energy LDOS spreads over 
the order of £o- Types of dispersion relation for the A- 
and B-phases are summarized in Table HI 

The spin degeneracy of the QPs is also different be- 
tween the A- and B-phases. In the A-phase, since the 
up-up and down-down spin Cooper pairs have the same 
chirality, the low energy QPs have the degenerate branch 
of the dispersion at the edge (Fig. GUJh)). I n the B- 
phase, since the up-up and down-down spin Cooper pairs 
have the opposite chirality, the low energy QPs have two 
branches of the dispersion at the edge (Fig. [Etc)). This 
difference manifests itself in the edge current, namely, 
the mass current in the A-phase (Fig. EJb)) and the spin 
current in the B-phase (Fig. [7]). 



B. Stripe phase 

We discuss the Majorana zero modes at a domain wall 
in the stripe phased We take the thickness of a film along 
the z-direction and the modulation of the OP along the 
x-direction. In the stripe phase, the fc z -component of the 
OP changes the sign at the domain wall perpendicular to 
the x-direction so that the pair breaking by the reflection 
at the surface of the film is prevented. Then, the OP is 

A^ky, Aj_fc z ) to the right 



described as A r j g ht = (A^k x 



of the domain wall and A 



left 



(Am**, 



-Aj_k z ) to 

the left of the domain wall, where A|| is finite everywhere 
and Aj_ vanishes at the domain wall. 

The QPs with the fc z -component across the domain 
wall have the finite energy Andreev bound states because 
they feel the sign change of the fc z -component of the pair 
potential, which is not the exact 7r-phase shift. In ad- 
dition, since the QPs with k z = feel the full gap of 
the pair potential, they are not excited in low energy. 
Therefore, the Majorana zero-energy QP is absent at the 
domain wall in the stripe phase. The domain wall is 
qualitatively different from the edge. 
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TABLE I. Types of dispersion relation at the upper or lower surface, side edge, and corner of a slab for 
the A- and B-phases. Dirac valley is formed from linear dispersion from E = along a certain momentum 
direction and dispersionless zero-energy modes along another momentum direction in long wavelength. 
Majorana cone is formed from isotropic liner dispersion from E — on the plane parallel to a surface in 
the momentum space in long wavelength, where the zero mode is in a point on the Fermi surface. The 
mark x indicates absence of zero-energy modes. 









upper or lower surface 


side edge 


corner 


A-phase 


thin slab (D 


= 860 


X 


Dirac valley 


Dirac valley 




thick slab (D 


= 


X 


Dirac valley 


Dirac valley 


B-phase 


thin slab (D ■ 


= 1460 


X 


Majorana cone 


Dirac valley 




thick slab (D 


= 30£o ) 


Majorana cone 


Majorana cone 


Dirac valley 



C. Experimental proposal 

There are several experimental means to detect the 
Majorana nature. Surface specific heat measurement, 
which was performed in connection with detection of the 
Andreev surface bound stated resolves the side edge 
contribution C sur f aco (T) = jT of the A-phase in a thin 
slab at low temperatures, where 7 oc N(E = 0), because 
the bulk contribution Cbuik oc T 3 which comes from point 
nodes where N(E) oc E 2 is distinguishable. Note that, 
in the A-phase, the contribution from the two upper and 
lower specular surfaces in the slab geometry is the same 
as that from the bulk. Thus, the surface specific heat 
C sm { acc (T) = jT of the Majorana QPs is distinctive. If 
the Z-vector direction is modulated near the side edge, the 
LDOS at the surface is N(E) = N(0)+N geip (E). The sur- 
face specific heat from the LDOS is also C sur f ace (T) = jT. 

In the B-phase, the zero-energy LDOS is localized at 
the corner of the order of 63. The contribution from the 
corner is interesting but smaller than that of the sur- 
face. Since we will discriminate Cbuik oc T~ 3 / 2 e~ A / ksT 
from the gap A and C sur f aco oc T 2 from the linear be- 
havior of LDOS N(E, r) oc E near E = 0, the Majorana 
fcrmion can be observed. In the slab with short thickness 
where the zero-energy modes are split at the upper and 
lower surfaces, the difference of the specific heat from the 
surface and bulk will not be distinctive because the gap 
structure and linear behavior of the LDOS coexist. 

The observation of the edge mass current in the A- 
phase, which is intimately connected with the intrinsic 
angular momentum^ is also hopeful. The magnitude 
of the edge mass current is unchanged in wider slabs 
than L = 40£o for which we have calculated. Consid- 
ering the observation of the torque from the edge mass 
current for a 10 mm x 7 mm x 0.6 fim slab sample 
(Bennett et al£§- had been used) by a typical torsional 
oscillator, the frequency shift from the edge mass current 
is of the order of 10 -23 KzM- The torque is too small to 
observe by a torsional oscillator because the magnitude 
is ~ NH, where N is the total number of 3 He atoms in 
the slab sample. We have to consider other experimen- 
tal methods. Also the edge spin current in the B-phase 
which flows thrce-dimensionally has been obtained quan- 



titatively. The techniques to detect the spin current is 
desired. The specific experimental proposal to observe 
the edge current is a future problem. 

The most direct evidence of the Majorana nature is 
derived from the observation of the anisotropic spin sus- 
ceptibility. If we use the Majorana nature of the edge 
state, in the A-phase, the local spin operators result in 
S x ~ S y ~ and only S z parallel to d-vector remains 
nontrivial for T <C T c *& This predicts the Ising-like spin 
dynamics for the local spin operator parallel to d-vector 
in the A-phase as well as that perpendicular to the edge 
in the B-phase ^ This is in sharp contrast to the suscep- 
tibility parallel to d-vector in the bulk A-phase which is 
suppressed at low temperatures according to the Yosida 
function. On the other hand, the susceptibility perpen- 
dicular to d- vector still assumes the bulk value, which 
is the same as it in the normal state. The anisotropic 
susceptibility has been discussed also by Shindou et alM- 
and has been calculated in the B-phase by Nagato et alZ 

QP scattering or QP beam experiments are extremely 
interesting. They were performed in the past on He 
where roton-roton scattering is treated^ and on the 3 He 
B-phase where the surface Andreev bound state is inves- 
tigated^^ Using this method, we may pick up Majo- 
rana QPs with a particular wave number. Particularly 
in the A-phase, the Majorana QPs from the edge is sep- 
arated from other QPs from the nodal region. 

Another option might be to use a free surface where the 
Majorana fermion surface state is formed. As shown by 
Kono^ it can be detected through the excitation modes 
of the floating Wigner lattice of electrons placed on the 
surface. We need a special, but feasible configuration of 
the experimental setups. 

Note that the recent work of transverse acoustic 
impedance measurements to detect the surface bound 
states in the superfiuid 3 He will derive the important 
information of the Majorana QPsi 49 ' 50 



VII. SUMMARY 

We have designed a concrete experimental setup to ob- 
serve the Majorana nature at the surface in the slab ge- 
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ometry. In connection with realistic slab samples, we 
have considered the upper and lower surfaces and the side 
edges including the corners with several thicknesses. We 
have demonstrated that the quasi-classical Eilenberger 
equation yields the quantitatively reliable information 
on physical quantities for the superfiuid 3 Hc A- and B- 
phases. Specifically, we have exhibited the difference of 
LDOS between the A- and B-phases and evaluated the 
mass current for the A-phase and the spin current for the 
B-phase quantitatively. Then, we have shown the influ- 
ence on the Majorana zero modes from the spatial varia- 
tion of /-vector for the A-phase in the thick slab and the 
energy splitting of the zero-energy modes for the B-phase 
confined in the thin slabs. The corner of the slab in the 
B-phase is accompanied by the unique zero-energy LDOS 
of corner modes. In addition, we have demonstrated the 
absence of the Majorana zero-energy QP at the domain 
wall in the stripe phase. On the basis of the quantitative 
consequences, it is proposed that the measurement of the 
specific heat, the edge current, and the anisotropic spin 
susceptibility provides feasible and verifiable experiments 
to check the Majorana nature. The control on the thick- 
ness of the slab is crucial to detect the Majorana surface 
states. The experiment controlling the thickness of the 
film of the superfiuid 3 He is interesting^ 
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where we describe the complex conjugate of the Matsub- 
ara frequency explicitly because that is important when 
we calculate LDOS. By the relations, we are allowed to 
sum only the positive ui n to calculate self-consistent pair 
potential and mass and spin currents. 

If we reverse the sign of relative momentum fc, the 
sign of spin-triplet components of the OP changes; on 
the other hand, that of a spin singlet component of the 
OP does not change. Specifically, general OP 



A(k,r) 



A tt (fe,r) A n (fc,r) 
A it (fc,r) A u {k,r) 



has a relation, A(— k, r) = — A T (fc,r), where a super- 
script T indicates transposition of a matrix. Since the 
sign of Fermi velocity also changes, v(—k) = —v(k), the 
quasi-classical Green's functions in particle-hole space 
have relations 



g(-k,r,u*) = g(k,r,u) n )*, 
g(-k,r,u)*) = g(k,r,u} n )*, 
/(-fc,r,0 = -/(fe,r,w n )*, 
f(-k,r,cj* n ) = -f(k,r,u ln )*. (A.2) 
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Appendix 

We use the symmetry of the quasi-classical Green's 
function in Eilenberger Eq. ([2]) to reduce computa- 
tional time. If we replace Matsubara frequency tu n with 
— w* , the quasi-classical Green's functions in particle- hole 
space have relations 

g(k,r,-w*) = -g{k,r,uj n )\ 
g(k,r, -w*) = -g{k,r,uj n )\ 
f(k,r,-u)*) = /(fc,r,w„) t , 
/(fc,r,-<) = /(fc,r, W „)t i (A.l) 



More reduction of computational time is possible by 
using mirror operators which define S x a = (—a x ,a y ,a z ), 
S z a = {a x ,a y ,~a z ), and S xz a = (-a x ,a y ,-a z ) with 
an arbitrary vector a = (a x , a y ,a z ). If the mirror opera- 
tors act on the center-of-mass coordinate of pair potential 
under the antiperiodic boundary condition mentioned in 
Sec. IIIIl A(fe, Sr) = A(Sk, r), where S is one among S x , 
S z , and S xz . Since Sv ■ SV = v ■ V, the quasi-classical 
Green's function satisfies 

g(Sk, Sr, w„) = g(k, r, uj n ). (A. 3) 



Therefore, the quasi-classical Green's function and pair 
potential are obtained self-consistently by numerical cal- 
culation for only the positive Matsubara frequency and 
one eighth of the Fermi surface in the coordinate —L/2 < 
x < L/2 and —D/2 < z < D/2 by the symmetry of the 
quasi-classical Green's function (|A. 1|) . (|A.2[) . and (|A.3[) . 
The calculation for the real space coordinates is carried 
out by parallel computing with OpcnMP. 
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